In situ measurement of the dynamic structure factor in ultracold quantum gases 
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We propose an experimental setup to efficiently measure the dynamic structure factor of ultracold 
quantum gases. Our method uses the interaction of the trapped atomic system with two different 
cavity modes, which are driven by external laser fields. By measuring the output fields of the 
cavity the dynamic structure factor of the atomic system can be determined. Contrary to previous 
approaches the atomic system is not destroyed during the measurement process. 
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I. INTRODUCTION 



The rapid experimental progress in the manipulation 
of ultracold quantum gases has enabled the creation of 
strongly correlated many-body systems, which are chal- 
lenging to describe theoretically [1]. The unambiguous 
identification of these phases requires precise measure- 
ments to characterize its properties. A powerful tool 
represents the response function, which provides static 
properties of the system as well as reveals information 
about the excitations of the examined state. Here, we 
propose a method for an in situ measurement of the re- 
sponse function to a external and weak probing field for 
a trapped gas of ultracold atoms. 

In seminal experimental measurements the Bogoliubov 
excitation spectrum in a Bose-Einstein condensate has 
been measured by probing the dynamic structure factor 
[2I, y]. Nowadays, this method has been applied to access 
the excitation spectrum within the BEC-BCS crossover 
[4, 5], as well as the smooth transition of the excitation 
gap from the superfiuid to the Mott insulating phase [6| . 
A major drawback of the present approaches is that the 
setup measures the number of excitations created in the 
systems: the atomic system is destroyed during the mea- 
surement process, and in addition, it requires the appli- 
cation of strong fields violating in most situations the 
condition of a weak probe. As a consequence, it becomes 
difficult to distinguish, whether the measurements probes 
the dynamical evolution of a strongly driven system or 
characterizes the ground state properties. Several alter- 
native methods to access information on ground state 
properties by measuring the light field passing through 
an atomic system have been discussed, e.g., the detec- 
tion of magnetic order |7H1Q| and particle fluctuations 
[ill, [l2|. Within a remarkable experiment, the possibility 
to access the dynamical structure factor by a measure- 
ment of the probe field has recently been demonstrated 
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FIG. 1: Proposed experimental setup: an ultracold atomic 
cloud interacts with two different cavity fields, which are 
driven by external lasers. The light fields leaving the cav- 
ity contain information on the dynamic structure factor of 
the atomic system, which can be recovered in a homodyne 
detection scheme. 



[l3| , which opens a way for weakly probing and analyzing 
cold atomic gases. 

In this article, we propose a setup for realizing an in- 
situ measurement of the linear response function for a 
cloud of trapped ultracold atomic atoms. The setup uses 
two different cavity fields, see Fig [TJ and is motived by 
recent progress of Bose-Einstein condensates coupled to 
high finesse cavities |l3-[l7|. Note that the experiment 
can also be performed in a ring cavity setup to avoid 
standing waves created by the probe lasers; similar se- 
tups have been proposed for cavity state preparation 
[iSl, [l9|. The main idea is that the interaction of the 
atomic system with the cavity modes transfers informa- 
tion on the response of the atomic system onto the cavity 
fields [sl, |Tl| . By measuring the fields that leave the cav- 
ities one can reconstruct the linear response function as 
well as the dynamic structure factor. Our method puts 



only modest requirements on the cavity finesse and al- 
lows for fast measurements with sufficiently high photon 
count rates without destroying the atomic system. 

For simplicity, we focus in the following on the density 
response function; the extension to alternative setups ac- 
cessing e.g., the spin structure factor is straight forward. 
For weak probe fields (/), the linear response function is 
defined by the deformation of the particle density 



(Pq(^)) = Xq(^)<?^q(^) 



(1) 



with the Fourier transformation relation (/)q(co') = 
J dtdr(p{t^r)exj[){iujt — iqr). While the real part of the 
response function x' accounts for the dispersive proper- 
ties of the media on the probe field, the imaginary part 
x" describes the creation of excitations within the media. 
Consequently, the response function contains important 
information about the two-particles excitations as well as 
collective excitations, which is obvious from its relation 
to the dynamical structure factor 



Xq(^) 



-7r[l — exp(— /3fia;)]5'(q, cj). 



(2) 



In previous experiments the dynamic structure factor 
5'(q, cj) has been studied extensively. However, it is im- 
portant to point out, that the real part and the imaginary 
part of the response function are not independent of each 
other, but rather are related via the Kramers-Kronig re- 
lation 



x;'H = — w/dw' 



,/2 _ ,,,2 



(3) 



As a consequence, it is possible to access information 
about the excitation spectrum by probing the real part 
of the response function off-resonantly and without cre- 
ating any excitations in the media. Such a method is in 
strong contrast to the current approach, where the cre- 
ated excitations strongly distort the probed state. The 
method present in this manuscript allows one to access 
both the real part of the response function as well as the 
imaginary part, and opens a way to study in situ and for 
very weak probes the response of an quantum many-body 
system with cold atomic and molecular gases. 



II. HAMILTONIAN DESCRIPTION 

We consider a system of N atoms with two internal 
states 1^) and |e) (see Fig. [2]), coupled to two cavity fields 
ttk and ttp with frequencies cjk and cjp, respectively. The 
dynamics of the cavity modes is governed by the Hamil- 
tonian 



He = huj\cal.a\i -\- Hujpa^^a^ 



(4) 



For a large detuning Ag of the cavity fields from the 
atomic resonance, we can adiabatically eliminate the 
state |e). Ignoring effects from spontaneous emission, the 
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FIG. 2: Internal level structure of the atomic system. Transi- 
tions between the ground state \g) and the excited state |e) are 
driven by cavity fields ^^k, ^p with a large detuning Ae from 
the atomic resonance. The momentum difference q = k — p 
of the cavity fields is transfered to the atomic system. 



interaction between the two cavity fields and the particle 
density operator pq of the atomic system takes the form 



Hi=g 



Po 



f aj^aic + ajjttj 



-q^k^p 



■ Pqfflpak 



, (5) 



with g being the two-photon coupling constant, while q = 
k — p denotes the difference in momenta between the cav- 
ity fields. The first term in Hj accounts for the Stark shift 
of the atoms and only leads to a renormalization of the 
cavity resonance frequencies, i.e., a;k,p = ^k,p + Poq/^- 
Note that here we assume po denotes the total number of 
atoms in the cavity and is a conserved quantity. In addi- 
tion we would also like to point out that the cavity lasers 
can also be focused and probe only parts of the atomic 
system. Then, po accounts only for the atoms within the 
mode volume of the lasers, and the Stark shift becomes 
a local shift in the chemical potential. 



III. WEAK COUPLING EXPANSION 

A. Free field solution 

We consider the terms involving the density operator 
p±q as small perturbations and solve the system first 
in absence of these coupling terms. Then, the quantum 
Langevin equation for the cavity fields reduce to 



dk = -ic^k<^k - -ttk + Vr^k (^) 



(6) 

(7) 



where ^kp(t) are the input fields of the cavities [20]. We 
assume that the cavity is driven by two lasers with fre- 
quencies r^k and r^p, i.e., (^kp(^)) = /5k,pexp(— ir^k,p^)- 
Then, the steady state solution for the averaged fields 
Q^k,p = (<^k,p) is given by 



ak(t) = 
ap{t) = 



^/7 



7/2 + [{tdk - ^k) 
7/2 + i(wp - i^p) 



/3kexp(-i!^kt) (8) 

/3pexp(-il^pi). (9) 



To simplify the analysis, we set the cavity detunings to 



ft. 



Qp = and assume equal 



zero, A = cjk 

driving fields \/3\^\ = |/3p|, while the generalization to un- 
equal driving fields is straightforward. Furthermore, we 
can absorb a phase factor in the definition of /3k, p such 
that a\^{t) = ao ex.p{—iQ]^t) and ap{t) = ao exp(— il^pt) 
withao = 2|/3k|/V7. 



B. Effects of fluctuations 

To study the effect of the atom-field interaction, we 
consider small fiuctuations around the free field solution. 
For this, we define a new set of operators as 



1 



A = ^ (akc^^^^ + ape 
B = ^ (a.e 



Anr,t\ 



iQi^t 



a^e 



iQr,t\ 



with the inverse transformation 



i^k^ 



ap = ^{A - B)e-'"p*. 



(10) 
(11) 

(12) 
(13) 



These operators have the expectation values (A) = 
2Y^2/7/3k and (B) = 0, and we expand them into their 
mean values and small fiuctuations SA^ 5B. Then, the 
interaction Hamiltonian Hj = Hi -\- H2 -\- Hs reduces to 



Hi 

H2 
H3 



/ dxp(x)l^ cos(qx — ujt) 



(14) 



iujt 



9eE 

2 



f/eff / 



(Pqe" 



-pV')(^^ + ^^0 (15) 
p-^e'^') {SB - SB^) (16) 



with a; = l^k — ^p, the potential strength Vq = ^|(A)p, 
and the effective coupling ^eff = ^1(^)1- The first term 
describes a classical driving field V^xt = Vo cos(qx — ut) 
for the atomic system, while the last terms account for 
the coupling between the atomic system and the cavity 
fields. For cj = 0, this set of equations reduces to the 
case previously studied in the context of cavity cooling 
|21| . We now introduce the quadrature operators 



Xb = 1 [SB + SB^ Pb 

Xa = \ [SA + Mt] Pa 



with the commutation relations 



-i [8B - SB"^] (17) 
-i [5A - 5A^] (18) 



[Xa,Pa] =[Xb,Pb]= i 
[Xb,Pa] =[Xa,Pb]= 0. 



(19) 
(20) 



The equation of motions for the cavity field then reduce 
to 



Xa = 
Pa = 
Xb = 
Pb = 



9eS 



[Pqe" 



Aujt] 



.^eff 



tfen r 



-iujt 



P-cie 



Pb- 



-^ 


(22) 


2x. 


(23) 




(24) 



These equations describe the back action of the atomic 
system with the particle density pq onto the cavity fields, 
i.e., a fiuctuation in the particle density operator of the 
atomic system Pq(t) infiuences the cavity fields within 
characteristic quantum- non-demolition setup. 



C. Linear response regime 

The leading term Hi drives a perturbation onto density 
of the atomic system with the field ^(t, x) = Vq cos(qx — 
ujt). The condition of a weak probe field reduces to l^ < 
Eg with £^5 ~ h^G^ /m the characteristic energy scale of 
the atomic system [22|. Then, the response of the atomic 
system to this external probe is well described within 
linear response theory 

W,,x))=/.<'.x'„*-*',x-xW.x'), (25) 

which implies for the present drive with momentum 
transfer q and frequency uj the response (pq(cj)) = 
VbXq('^)/2. Then, it is possible to replace the operators 
in the Langevin equations by their expectation values, 
and we obtain the coupled differential equations 



i^XA) = -'-(Xa) 

±iPA) = -^-fVox'^iu:) -1{Pa) 



di 
^AXb) = -^-^V,^:^{u:)-1{Xb) 



h 

geff 

2h 



At 



dt'~"' 2 

We can immediately identify the steady state solution 



(26) 

(27) 
(28) 
(29) 



(Xa) 
{Pa) 

(Xb) 



{Pb) = 

2geffVo 

h-f 

gesVo 



xM 



hj 



x"(^). 



(30) 
(31) 

(32) 



Consequently, the linear response of the atomic sys- 
tem is imprinted onto the cavity fields: the imaginary 
part of the response function describing the creation of 
two-particle excitations and collective excitations in the 
atomic system is encoded by a shift in the amplitude op- 
erator Xb accounting for the scattered of photons from 



one cavity mode onto the other. In turn, the real part of 
the response function characterizing the dispersive part 
of the media leads to a shift in the phase Pa of the fields. 



IV. EXPERIMENTAL PARAMETERS 

Finally, we provide the experimental parameters for 
the measurement of the response function. We will focus 
on the measurement of the phase quadrature (Pa), as it 
exhibits two advantages over the measurement of (Xb)'- 
First, we note that (Pa) is directly given by the sum of 
the expectation values (Pk) and (Pp) of the physical cav- 
ity fields ttk and ap, which can be efficiently measured 
within a homodyne detection. For appropriate phases 
of the cavity fields this means that the measurements 
are done against zero background, which is not the case 
when measuring (Xb)- Second, the Kramers-Kronig re- 
lation Eq. (J3]) allows us to compute Xq{^) on resonance 
by using measurement data taken off resonance. There- 
fore, we are able to obtain the dynamic structure factor 
of the atomic system without transferring energy onto 
the atomic system. 

The validity of the expansion and ignoring the 
quadratic terms in the fiuctuating fields requires a large 
photon number in the cavity, i.e., {A) ^ 1. In turn, 
the weak driving field Vq ^ Eg implies a weak coupling 
g = Vo/\{A)\'^ < Es/\{A)\'^. Such weak coupling is eas- 
ily achieved by controlling the detuning to the excited 
level, see Fig. [2l Note, that his weak coupling is in high 
contrast to the recent experimental works on high finess 
cavities [14h17i]. and significantly simplifies the require- 
ments for the experimental realization of such a setup. 
The phase quadratures can be efficiently been measured 
within a homodyne detection with the sensitivity, see 
Fig.S 



S- 



{Pa) (Pk> + (Pp) 



(^> 



{^> 



(33) 
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FIG. 3: Phaser diagram illustrating the measurement of (Pa) 
as the sum of the physical observables (Pk), (Pp)- A signal 
can be detected only when the ratio {Pa)/{{X\^) + (Xp) is 
large compared to phase fluctuations. 



The range of interesting values for the response function 
are Xq(^) ^ X/Es with X the total number of parti- 
cles and Eg the characteristic energy scale of the atomic 
system. This requires a sensitivity to measure the phase 
angle 



X^V^ 

n-i Es 



Vo X Fo 

h^\{A)\^Es 



(34) 



For the last equation, we have used the relation between 
the driving ffeld Vq = g\{A)\'^ and the number of photons 
in the cavity. The number of atoms within the cavity 
are typically in the range X ~ 10^, and therefore, with a 
cavity decay rate 7 ^ 100 GHz corresponding to a qual- 
ity factor Q ~ 10^ with up to 100 photons within the 
cavity, and a characteristic energy Eg ^ 1 KHz, the re- 
quired sensitivity is in the range S ~ 10~^, which can be 
achieved with present techniques [23|. 



V. CONCLUSIONS 

In summary we have demonstrated the experimental 
feasibility of an in situ measurement of the dynamic 
structure factor of arbitrary atomic systems, providing 
a novel tool for the coherent manipulation of ultracold 
quantum gases. The proposed method employs two cav- 
ity ffelds which are driven by external lasers. We have 
shown that the required resolution can be achieved using 
a cavity finesse that has already been reached in present 
experiments. Furthermore, the method allows for suffi- 
ciently fast measurements. 

The current setup is in close analogy to the recently 
experimentally realized system [13]. The main difference 
is, that in the experimental setup a strong imbalance 
between the two photon modes is applied. As a conse- 
quence, the validity of the expansion up to linear order in 
the ffuctuating fields is satisfied by the strong coupling 
laser, while for the weak probe beam a very low pho- 
ton number can be used. Then, the present calculations 
can be performed again in a straightforward manner with 
the only modification, that the real part of the response 
function is written on the phase quadrature Pp , while the 
imaginary part is encoded onto the amplitude quadrature 
Xp. In the experiment the amplitude quadrature was 
accessed with a heterodyne detection [13], while the real 
part of the response function again could be probed by a 
simpler homodyne detection. We would like to point out 
that the contrast of the signal can be strongly enhanced 
by using a cavity for the weak probe field. 

Finally, we would like to point out, that the present 
setup can also be used to study the relaxation of exci- 
tations in the atomic system: first, the lasers imprint 
several excitations on resonance into the atomic system, 
while at a later stage the probe of the real part allows 
one to analyze whether these excitations are still present 
or have relaxed. 
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